Microwave Spectroscopy of a Cooper-Pair Transistor Coupled to a Lumped-Element Resonator 
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We have studied the microwave response of a single Cooper-pair transistor (CPT) coupled to a lumped- 
element microwave resonator. The resonance frequency of this circuit, f r , was measured as a function of the 
charge n g induced on the CPT island by the gate electrode, and the phase difference across the CPT, which 
was controlled by the magnetic flux in the superconducting loop containing the CPT. The observed //•("[>, <fe) 
dependences reflect the variations of the CPT Josephson inductance with n g and <j>B as well as the CPT exci- 
tation when the microwaves induce transitions between different quantum states of the CPT. The results are in 
I excellent agreement with our simulations based on the numerical diagonalization of the circuit Hamiltonian. 

This agreement over the whole range of n g and <ps is unexpected, because the relevant energies vary widely, 
from O.IK to 3K. The observed strong dependence f r (ng,<j>B) near the resonance excitation of the CPT provides 
a tool for sensitive charge measurements. 
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The Cooper-pair transistor (CPT) is a three-terminal device which consists of a mesoscopic superconducting island connected 
to two leads by two Josephson tunnel junctions (JJs) (see, e.g. U EJ and references therein). The behavior of this device is 
controlled by two energies: the charging energy per junction, £'c=e 2 /2Cy ( Cj is the capacitance of a single tunnel junction), and 
the Josephson coupling energy Ej. The energies Eq and Ej could be made of the same order of magnitude by reducing the tunnel 
O junction in-plane dimensions (typically, down to 100 — 200nm for Al —AlO x —Al junctions). The energies of quantum states of 
^ the CPT are 2e periodic in a continuous charge n g = C g V g /e induced on the island by a capacitively coupled gate electrode. Here 
Q-i C g is the capacitance of the capacitor formed by the island and the gate electrode, V g is the voltage applied to this capacitor. The 
sensitivity of the CPT characteristics to the induced charge makes this device a very sensitive electrometer which, in particular, 
can operate in a low-dissipation dispersive mode (3]|4|. The interplay between the Josephson effect and Coulomb blockade leads 
to a quantum superposition of charge states in the CPT, which forms the basis for quantum computing with superconducting 
charge qubits [5-7 1. Since the first demonstration of the coherent superposition of states in the CPT more than a decade ago, the 
I CPT has been used as a test bed for many novel experimental techniques employed in the research on superconducting qubits. 

The microwave experiments with CPTs can be broken down into two main categories. In the first type of measurements, the 
CPT remains in its ground state because of a large mismatch between the probe signal frequency and the excitation frequencies 
of the CPT. During this adiabatic operation, the CPT can be described by its effective microwave impedance. This impedance, 
depending on the parameters of the Josephson junctions and the coupling of the CPT to the readout circuit, could be predom- 
inantly inductive (the Josephson inductance, the second derivative of the CPT energy in phase [8]) or capacitive (the quantum 
capacitance, the second derivative of the CPT energy in charge ll9l- fTTIl ). Note that if the CPT is coupled to a resonator and 
their levels are close in energy, the entanglement of the CPT and resonator states affects the impedance of this circuit even if 
the microwaves do not induce transitions between the states. In the latter case, the impedance-based description of the CPT is 
insufficient, and the solution of the quantum Hamiltonian of the system «CPT + read-out circuit* is required. In the second type 
of measurements, the microwaves induce transitions between different quantum states of the CPT. This, in particular, enables 
the preparation and manipulation of coherent superpositions of the ground and excited states in the quantum-computing-related 
applications of the CPT. 

In this paper, we present the microwave spectroscopic study of a CPT which probes both the ground state and excited states 
y—i of the CPT over wide ranges of the charge n g and the phase difference across the CPT, The phase was controlled by the 
magnetic flux in the superconducting loop containing the CPT. The CPT microwave response was analyzed by measuring the 
• i-H resonance frequency /,- of a lumped-element microwave resonator coupled to a CPT. A relatively high quality factor of this 
circuit allowed us to perform the measurements in the low-power regime with an average number of photons in the resonator 
less than one. When the detuning between the microwave frequency and the excitation frequencies of the CPT was large, the 
dependence /,•(«£, <fe) mostly reflected the variations of the CPT Josephson inductance with n g and <j>B- On the other hand, an 
avoided crossing of the CPT and resonator levels was clearly observed when the CPT excitation frequency was tuned to the 
resonator frequency by varying n g and 0g. The strong n ? -dependence of the circuit response in this regime provides a tool for 
sensitive charge measurements. The overall dependence f r (n g ,^B) is in excellent agreement with the simulations based on the 
numerical diagonalization of the circuit Hamiltonian. This agreement provides a stepping stone for the understanding of more 
complicated superconducting circuits intended for quantum computing, including multi-junction circuits envisioned as protected 
superconducting qubits fT2l \l3l . 

The paper is organized as follows. In Section [II] we describe the samples and measurement techniques. The details of 



numerical simulations of this circuit are provided in Section III The experimental results are discussed and compared with 
numerical simulations in Section iPVl 
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Figure 1: Schematics of the "CPT+LC resonator" circuit, (a) The on-chip circuit layout. The superconducting shunting wire ("shunt") serves 
two purposes: it reduces the coupling of the CPT to the LC resonator (and, thus, reduces external noises), and it forms, in combination with 
the CPT, a superconducting loop. The magnetic flux in the loop controls the phase difference across the CPT. (b) The circuit diagram used for 
modeling in Section 3. 



II. DEVICE FABRICATION AND MEASURING TECHNIQUES. 
A. Circuit Design 

The schematics of the tested circuit is shown in Fig{T] The CPT is inductively coupled via a narrow Al wire ("shunt") to a 
lumped-element LC resonator. The LC resonator consists of a meandered 2-jj.m-wide Al wire with L = 5 nH and an interdigitated 
capacitor (2-|im-wide fingers with 2 y,m spacing between them) with C = 100 fF. It was designed to be strongly coupled to the 
microstrip line in order to increase the signal-to-noise ratio in the low-microwave-power measurements with the average photon 
occupancy of the resonator n p h < 1. The typical values of the internal and loaded quality factors for these resonators (not 
coupled to the CPT) were 50,000 and 20,000, respectively. High Q values enable sensitive measurements of small changes in 
the microwave impedance of the tested device induced by the variations of n g and (jig. Outside of its bandwidth, the resonator 
efficiently decouples the CPT from external noises. An additional protection of the CPT from external noises is provided by 
the shunt: the kinetic inductance of this superconducting wire, L s i wnt — 0.5 nH, is more than 10 times smaller than the effective 
Josephson inductance of the CPT, which significantly reduces the phase fluctuations across the CPT. The LC resonator was 
inductively coupled to a 2-port Al microstrip line with a 50 f2 wave impedance. The gate electrode of the CPT is coupled to the 
central island of the CPT through a capacitor C g = 0.2 fF. 

The sample was mounted inside an r/-tight copper box which provided the ground plane for the microstrip line and LC 
resonator. This box was housed inside another r/-tight copper box in order to attenuate stray photons which may originate from 
warmer stages of the cryostat [ 14 1. This nested-box construction was anchored to the mixing chamber of a cryogen-free dilution 
refrigerator at a base temperature of 20 mK. 

B. Device Fabrication 

The Cooper-pair transistor, the lumped-element LC resonator, and the microstrip line were fabricated within the same vacuum 
cycle using multi-angle electron-beam deposition of Al films through a nanoscale lift-off mask. To minimize the spread of the 
junction parameters, we have adopted the so-called "Manhattan-pattern" bi-layer lift-off mask formed by a 400-nm-thick e-beam 
resist (the top layer) and 50-nm-thick copolymer (the bottom layer) (see, e.g., IT5l[T6l and references therein). In this technique, 
Josephson junctions are formed between the aluminum strips of a well-controlled width overlapping at a right angle. After 
depositing the photoresist on an undoped Si substrate and its patterning with e-beam lithography, the sample was placed in an 
ozone asher to remove any traces of the photoresist residue. This step is crucial for reducing the spread injunction parameters. 
The substrate was then placed in an oil-free high-vacuum chamber with a base pressure of5x 10~ 9 Torr. The axis of the rotatable 
substrate holder forms an angle of 45° with the direction of e-gun deposition of Al. During the first Al deposition, the substrate 
holder was positioned such that only the "avenues" were covered with metal; no metal was deposited in the "streets" because the 
mask thickness is greater than the width of the "streets". Since the mask thickness is 0.45jUra, this technique is suitable for the 
fabrication of JJs with lateral dimensions up to 0.3 x 0.3/Xm 2 . The thickness of this first Al film, which forms the central island 
of the CPT, was 20 nm. Without breaking vacuum, the surface of the bottom electrodes was oxidized at 100 mTorr of dry oxygen 
for 5 minutes. After evacuating oxygen, the substrate holder was rotated by 90°, and 60-nm-thick top electrodes were deposited 
along the "streets" (no aluminum is deposited in the "avenues" at this stage). The central island of the Cooper-pair transistor 
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was always deposited during the first Al deposition, and its thickness was smaller than that of the leads; this is important for 
preventing the quasiparticle poisoning [17]. Finally, the sample was removed from the vacuum chamber and the lift-off mask 
was dissolved in the resist remover. The spread of the resistances for the nominally identical JJs with an area of 0. 15 x 0. 15/im 2 
did not exceed 10%. More than 10 devices with Ej/Ec = 1 .5 — 3 have been studied and the results have been successfully fitted 
with the numerical simulations; below we discuss two representative samples. 

C. Measurement Technique 

The microwave response of the coupled system " CPT + LC resonator" was probed by measuring the amplitude and phase 
of the microwaves traveling along a microstrip line coupled to the resonator. Figure [2] shows a simplified schematic of the 
microwave circuit, which is similar to the one used in work [ 18 1. The cold attenuators and low-pass filters in the input microwave 
line prevented leakage of thermal radiation into the resonator. On the output line, a combination of low-pass filters and two 
cryogenic Pamtech isolators (-18 dB isolation between 3 and 12 GHz) anchored to the mixing chamber were used to attenuate 
the 5 K noise from the cryogenic amplifier. The DC line for the gate voltage control was heavily filtered with a combination of 
room temperature LC and low temperature RC filters, followed by a stainless steel powder filter, and a 1:1000 voltage divider. 

The probe signal at frequency 0)2, generated by a microwave synthesizer (Anritsu MG3694B), was coupled to the cryostat 
input line through a 16 dB coupler. This signal, after passing the sample, was amplified by a cryogenic HEMT amplifier 
(Caltech CITCRYO 1-12, 35 dB gain between 1 and 12 GHz) and two 30 dB room-temperature amplifiers. The amplified 
signal was mixed by mixer Ml with the local oscillator signal at frequency (0\, generated by another synthesizer (Gigatronics 
910). The intermediate-frequency signal a(t) — asin(J7f + <p) +noise(t) at £1 = (coi — «2)/2ti = 30 MHz was digitized by 
a 1 GS/s digitizing card (AlazarTech ATS9870). The signal was digitally multiplied by sin(Slf) and cos(Slf), averaged over 
many (typically, 10 6 ) periods, and its amplitude a (proportional to the microwave amplitude S21) and phase <p were extracted 

as a = J (a 2 (/)sin 2 fi? + a 2 (f)cos 2 Of) and <p = arctan ((a 2 (f) sin 2 £2f) / (a 2 (f)cos 2 £2f)), respectively (here (...) stands for the 
time averaging over integer number of periods). The value of <p randomly changes when both (0\ and (O2 are varied. To eliminate 
these random variations, we have also measured the phase <po of the reference signal provided by mixer M2 and digitized 
by the second channel of the ADC. The phase difference <p — <po, being dependent at fixed n g and 0# only on the electric 
length difference between the microwave lines inside and outside of the cryostat, is immune to the phase jitter between the two 
synthesizers. The low noise of this setup allowed us to perform measurements at microwave excitation level down to -140 dBm 
which corresponded to sub-single-photon population of the tank circuit. 

III. MODEL HAMILTONIAN AND NUMERICAL SIMULATIONS 

We begin with the discussion of the theoretical model of a more general circuit which contains an arbitrary Josephson device 
coupled via a superconducting "shunt" to a microwave resonator. The only limitation on the device parameters is that all 
characteristic energies of the device are much smaller than the effective inductive energies of all superconducting wires in the 
circuit, Et i = ti 2 /(2e) 2 Lj. The resonance frequency of the circuit might be of the same order or even very close to the device 
excitation energies, which would lead to the level repulsion. To simplify the notations, we shall use below the units h = 2e = 1 
(e.g., in these units V = d(j) jdt) and restore the physical units at the end, where we apply this model to the specific case of a 
device that consists of two Josephson junctions and one superconducting island, i.e. the CPT. 

The generalized circuit shown in Fig. lb includes two loops: the long meandering wire, the shunt, and a large capacitor C form 
one loop (referred below as the "resonator" loop), and the device and the shunt form another loop (referred as the "device" loop). 
The inductance of the meander (shunt), and the phase difference across this element are L m (L s f,) and § m (0 S /,), respectively. The 
difference between the device phase <po and the shunt phase (j)^ is due to the time-independent magnetic flux <t> in the device 
loop: 0o — §sh = (t>B, where <ps = 27£ < t , / < t > o, 't'o is the flux quantum. The voltage drops across the device and the shunt are equal: 
Vsh = Vd- The whole circuit is described by the Lagrangian 

~S? = T sh (V sh )+T m (V m ) + ^(V sh +V m ) 2 -^E sh (l) 2 h -^E m <j) m 2 (1) 
+ J%(0o,Vo). 

Here T m (V) (T s h(V)) is the time-dependent part (i.e. the kinetic energy) of the response of the meander (shunt), and Jz?£)(0o, Vq) 
is the device Lagrangian which also depends on the internal degrees of freedom (phases) of the device. In the BCS theory, 
the response of a superconducting wire with the static energy El has a scale of El at all frequencies 03 < A; it is a function 
of the dimensionless parameter V/A: T = El/(V /A) = (1 /16)(V/A) 2 + 0((V /A) 4 ) where A is the superconducting gap. This 
equation implies that at low frequencies, the wire impedance acquires, in addition to the kinetic inductance, a small capacitive 
component C' = El/ (8A 2 ) (these capacitances are shown in Fig. lb by dashed lines). 
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Figure 2: Simplified circuit diagram of the measurement setup. The microwave signal at o>2 transmitted through the microstrip line coupled to 
the "LC resonator + CPT" circuit is amplified, mixed down to an intermediate frequency ©1 — 0)2, and digitized by a fast digitizer (ADC). The 
second channel of the ADC is used to digitize the signal from an additional mixer (M2), which provided the reference phase <po (see the text). 
The gate voltage V g is applied to the capacitor C g using a heavily filtered DC line. 



We shall assume that the device Lagrangian is given by the sum of the Josephson and electrostatic energies: 

&d(<I>d,V d ) = l^QjViVj + 'EJij cosifr-Q j-Qij) (2) 

Here phases 0, and corresponding potentials V, describe both the internal degrees of freedom of the device and the shunt phase. 

Because the potential energy of the resonator loop is much greater than that of the device, the effect of the device on this loop 
can be treated as a small perturbation. In the abse nce of the device, the resonator loop has two modes: the harmonic oscillation 
of the total phase S /, + <j) m with frequency ©0 = yl/ L/C where L' = E m l +E~j}, E^sh = E m (j) m , and the orthogonal mode with 
<p m + (ps/j « 0. Because the large capacitance C does not participate in the second mode, the frequency of this mode is determined 
by the superconducting gap, the only energy scale in this case. We shall assume that A 3> (Oq so both real and virtual excitations 
of this mode can be neglected as well as the contribution of the capacitances C' sh and C' m to the effective capacitance C of the first 
mode. Note, however, that the virtual processes involving the second mode are small only in COq/A and might not be completely 
negligible in a realistic situation. If these processes are neglected, there is only one relevant degree of freedom, the phase across 
the device 0o- The effective Lagrangian is reduced to 

■&eff = - - 0b) 2 + ^d(^V Q ) (3) 



with C L = C(l +E sh /E m ) 2 and 
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E L = E sh {\+E sh /E m ). (4) 
Fluctuations of the phase fa in the low-energy states of the oscillator mode are very small: 



((fa -fa) 2 ) = A 2 (2« + l)«l 
cop 
2E L 



A 2 



where n is the quantum number of the oscillator states. This allows one to replace the solution of the full problem by the solution 
of the simplified model in which we expand the interaction term in small phase fluctuations. 

It will be more convenient to use the Hamiltonian formalism in which the conjugated degrees of freedom are phases and 
charges. The total Hamiltonian is the sum of three parts, the Hamiltonians of the resonator (Hr), device (Hp), and interaction 
between them (i/,„ r ) : 

On n 1 ~ 

Hr = ^ql + -E L {fa~fa) 2 (5) 
H int = C L l £ qoCojCJ^qji-nfi-YtJiocosifo-Qj-c&oi) (6) 

i,j>0 i 

Hd = \ £ (qi - ni)Cr/ (q; - nj ) - - £ J i} cos(^ - fa - <t>,j) (7) 

L ij>0 L ij>Q 

The coupling to the inductor charge fluctuations contains the inverse of the capacitance matrix (~ C£ l ) and thus is very small. 
Thus, even though the charge fluctuations across the shunt are not small: 

their effect on the coupling can be treated perturbatively. In the leading order in the interaction, we need to keep only two types 
of terms. The first type is quadratic in phase fa and diagonal in the basis of resonator states. The second type is linear in fa and 
off diagonal in this basis. The quadratic terms in the inductor charge are absent, so the charge coupling appears only due to the 
off diagonal terms that are linear in qo. The Hamiltonian equivalent to Q becomes 

He/f = {(Oo +2A 2 S)(a + fl + 1/2) + (AJ+ —Q) a + h.c. (8) 

Here of (a) is the creation (annihilation) operator for the harmonic oscillations of the circuit, J, Q and S are operators acting on 
the device which forms are obtained by expanding the interaction Hamiltonian 



dL 
d<j>o 



= -JViosin(0j+<I>o,-) (9) 

o i 

Q =C L l £ CojCJ^qji-rii) (10) 

1 " L = ^y i0 cos(^+4> ,) (ii) 

=0 1 i 



2 d(j) ( 



We now estimate the scale of the frequency deviations induced by these perturbations. In th e natura l units of the resonator 
frequency a>o, the scales of the perturbing operators are AJ/cOq ~ Ej / yj (QqEl, A _1 Q/oo ~ \J (Oq/El and A 2 Z,/cOq ~ Ej/El- 
The operator E is diagonal in the oscillator states, so it directly results in the frequency shift 8(Qz/(Qq ~ Ej/El. Because the 
non-diagonal elements affect the level of the resonator only in the second order of the perturbation theory, the effect of the J and 
Q operators depends on the gap between the levels in the combined resonator/device circuit. Far away from the full frustration 
and charge degeneracy point (fa = K, n g = 0.5), the device is characterized by large Ej^> coq and the energy levels are separated 
by large gaps, so the smallest gap is due to the resonator: 8E = (Oq. In this case the frequency shifts are 8(Oj/(Oq ~ E 2 /{(OqEl) 
and 8coq/cOq ~ COo/Ei respectively, which implies 8(Qj ^> 8(Oz ^ 8cOq. The effect induced by the phase and charge coupling 
grows when the gap between the levels coupled by these operators becomes small, but the phase coupling remains larger than 
the charge coupling for the devices with Ej (Oo- This increase of the frequency shift occurs, for instance, when the device 
level crosses the first resonator level. 
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Figure 3: (color online) Panel a shows the frequencies of three lowest-energy levels for the coupled system "CPT + LC resonator" plotted as a 
function of n s at a constant phase <j>s = 0.97k. The model parameters are Ej = hx 17.2GHz, Eq = hx 8.6GHz and El — hx 5650GHz. In 
this very frustrated regime, the lowest level of the CPT crosses the lowest level of the resonator with approaching the charge degeneracy point 
n g = 0.5, which results in the avoided level crossing. Panel b shows the blow-up of the theoretical curves in the region of avoided level crossing 
(black curves), together with the dependence of the color-coded microwave amplitude S21 on the microwave frequency / and n g measured for 
one of the tested CPTs at <ps = 0.97k. The minimum of the microwave amplitude corresponds to the resonance frequency. The dependence of 
S21 (/) at the charge degeneracy point n g = 0.5 is shown in Panel c. 



We now write down the explicit equations for the Cooper pair box. In this case, the internal degrees of freedom are limited to 
one phase (and the conjugated charge). Assuming equal capacitances and Josephson energies of the CPT junctions, we have 



H C pt = 4£ c .( ?1 -« s ,) 2 -^[cos(0i)+cos(^ 1 +^)] (12) 

J = -Ejamifa+fo) (13) 
2C 

Q = —E c { qi -n g ) (14) 

E = -Ejcos(<j)i+(j) B ) (15) 

Here we restored the physical energy units E c = e 2 /2Cj. For practical computations it is sufficient to retain the first few levels 
of the resonator (a' a < n mux = 3) and some number, ng, of the charging states. The Hamiltonian ^ becomes 3ng x 3ng matrix. 
Because the wave function of the charge decreases exponentially at large charges, ^(q) ~ exp(— yjE c /Ejq 2 ), it is sufficient to 
consider ng ~ 10 for accurate computations. The straightforward numerical digonalization of the Hamiltonian ([Sji leads to the 
theoretical predictions that can be compared with the data. 

Our experimental situation corresponds to E c ~ 2fi©o and Ej ~ 4nt0o. In the absence of frustrations, the frequency of the 
lowest CPT level is very high: (O p = \/%E c Ej ~ lOtOo- The frequency of the lowest CPT level decreases as the magnetic field 
frustrates the Josephson coupling and/or with approaching the charge degeneracy (n g — 0.5). Figure [3^ shows three low-energy 
levels of the system "CPT+resonator" with the parameters typical for our experiment. Note that for the studied circuits, only 
the combined effect of flux- and charge-induced frustrations brings the frequency of the first device level below that of the 
resonator, otherwise the device resonance frequency significantly exceeds that of the resonator even at full flux frustration (e.g. 
(O = 2E C > 4fflo at n g = 0). 



IV. EXPERIMENTAL RESULTS AND DISCUSSION 



Below we present the measurements of the amplitude S21 of the transmitted microwaves (unless otherwise specified) at the 
base temperature T = 20 mK. Most of the data (with the exception of the data in Fig.3b,c) are shown for only one representative 
device. The resonant dependence of S21 on the microwave frequency / , measured for this device at 0b = and n g = 0, is shown in 
the inset to Fig. [4] The resonance frequency depends periodically on n g and 0b; for example, the dependence f r (n g , <fe) measured 
at n g = is shown in Figure 0] The period in charge is An g — 2e at the base temperature (see Fig. 5); it changes from 2e to e at 
higher temperatures (>300 mK) due to the presence of thermally excited quasiparticles. Note that the total time of acquisition 
for the data shown in Fig. 3b was approximately 20 minutes; over longer time intervals, the periodicity of fr(n g ,$B) might be 
disrupted by the motion of non-equilibrium quasiparticles to/from the CPT island (the so-called "quasiparticle poisoning") ifTTll 
or other types of charge fluctuations fl9l . The high stability of the charge on the CPT island indicates that (a) the combination 
of a larger superconducting gap of the CPT island and its relatively large charging energy protects the CPT from quasiparticle 
poisoning, and (b) the double-wall //-tight sample box shields the device from stray high-energy photons. The microwave 
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Figure 4: Dependence of the resonance frequency f r (n g = 0.<Pb) on the magnetic flux <S, which controls the phase difference across the 
CPT, <j>B = 2^<J>/<I>o. The solid curve shows the numerical simulation with the fitting parameters discussed in the text. The inset shows the 
dependence of the microwave amplitude S21 on the frequency near the resonance at <j>g = and n g = 0. 

photon energy E p i, w/iX7 GHz is insufficiently large to excite the CPT at n g = 0: indeed, according to our simulations, the 
lowest excitation frequency for this device exceeds 30 GHz even at full flux frustration (O/Oo = 0.5). In this case the variation 
of the resonance frequency f r with magnetic flux reflects the ^-dependence of the CPT impedance in its ground state. 

The dependences of the resonance frequency on n g are illustrated by Figs. |5la-e), where the color-coded microwave amplitude 
52i is plotted versus / and n g for several values of the magnetic flux in the device loop. The black curves in Fig. [5] show 
the results of fitting the experimental data with our numerical simulations. All these curves were generated with the same 
set of fitting parameters: Ec = h x 16 GHz, Ej = h x 32 GHz and Ei = hx 5720 GHz (note that not only the amplitude of 
the resonance frequency modulation, but also the absolute values of f r are pre-determined by these parameters). The fitting 
procedure is very sensitive to the choice of these parameters: we believe that they are determined with an accuracy better than 
10%. The extracted charging energy coincides (within 5% accuracy) with an estimate of Ec based on the junction area, the 
specific geometrical capacitance for Al tunnel junctions (50 fF//im 2 , see e.g,. |20|), and the electronic capacitance of Josephson 
junctions, C e = 3/16(/?g//?)e 2 /A (0.3 fF at R — 3 kil) 12T1I221 . The Josephson energy estimated on the basis of the Ambegaokar- 
Baratoff relationship [23 1 using the normal-state resistance of a test junction deposited on the same chip is -40% greater than 
the fit value of Ej. 

Generally, one expects that Josephson circuits can be accurately described by the Hamiltonian consisting of Josephson and 
charging energies (cf. Eqj2]i only if all energy scales are smaller than A. Away from full frustration, the energy of the CPT 
excited state is of the order of Josephson plasma frequency ps 3.2K, which is comparable to A. Thus, the excellent agreement 
between the experimental data and numerical modeling, observed over the whole range of n g and 0b, is quite surprising. 

It is worth noting that the circuit modeling based on the numerical diagonalization of the circuit Hamiltonian is essential for 
fitting the data for devices with Ej/Eq ~ 2. For example, the analytical solution for the Josephson inductance in the CPT ground 
state, calculated within the two-level approximation (Eq. 4 in Ref. J8j), overestimates the amplitude of the /,- (n g ) dependence at 
small flux frustrations by almost an order of magnitude. The latter solution provides more accurate fitting of the experiment at 
larger frustrations (<pg ~ (0.8 — 9)n), but becomes inadequate again at 0# > 0.9471 when the avoided level crossing is observed. 

The evolution of these dependences reflects the modification of the CPT spectrum with n g and <ps ■ For a small phase difference 
(j)B (Figs|5^-b), the lowest CPT excitation frequency well exceeds the microwave frequency (which is close to the resonance 
frequency of the resonator), and the CPT remains in its ground state for all n g including the charge degeneracy point (n g = 0.5). 
In this regime, the dependences S2t{f,n g ) mostly reflect the variations of the CPT impedance with n g in the CPT ground state. 
At a larger frustration §b = 0.9n (Fig. |5J;), the lowest CPT level approaches the lowest resonator level n g = 0.5 (but has not 
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crossed it yet), and the entanglement of the qubit and resonator states leads to a complicated shape of the f r (n g ) dependence. 
Finally, with the further approach to full frustration (0g > 0.947T, Figs. [5tl,e), the CPT excitation frequency becomes smaller 
than the resonance frequency of the LC resonator, and the shape of the fr[n g ) dependences abruptly changes: they are strongly 
affected by the avoided level crossing. The S%i (f>n g ) plots in this regime consist of two sets of curves. The lower set of curves 
corresponds to the lowest energy level of the combined system "CPT + resonator" (this level coincides with the CPT lowest level 
when approaching the charge degeneracy points, i.e. far away from the resonance frequency of the resonator). The upper set of 
curves corresponds to the first excited level of the system "CPT + resonator": with approaching the charge degeneracy point, 
this level descends from higher energies to its lowest position at n g = 0.5. The visibility of the upper set of curves depends on 
the proximity between the CPT and resonator levels. Indeed, if the energy of the CPT resonance at n g = 0.5 is much lower than 
the first resonator level, the upper-curve "cone" is very sharp, and the corresponding microwave resonance is smeared even for 
small deviations of n g from 0.5 (this case is illustrated by Figs. |5jl,e). On the other hand, the "cone" becomes broader when the 
intersecting CPT and resonator levels are close to one another: in this case illustrated by Fig. lb, we were able to follow the 
upper set of curves over the frequency range of -15 MHz. 



For both devices, whose dependences S21 (/,«#) are shown in Figs. 3 andB] we observed a double-resonance structure at full 
frustration and charge degeneracy, depicted in Figs. 3b,c and 5e. The second(low-frequency) resonance appears as a "shadow" 
of the resonance observed at n g = 0. The appearance of this resonance, much weaker than that at n g = 0, implies that there 
are fluctuations of the island offset charge ±e which are fast at the measuring time scale -0.1 s. These fluctuations change 
the effective n g from 0.5 to 0. We attribute these fluctuations to the non-equilibrium quasiparticles moving between the CPT 
island and the leads. At n g = 0.5, the energy of a quasiparticle on the island exceeds the energy of quasiparticles in the leads 
by 8 A — (\/2)Ec- Here 8 A is the difference between superconducting gaps in the island and the leads due to the difference in 
the thicknesses of these Al films; we estimate 8A to be -kg x 0.3K. In our devices, the probability of these fluctuations is small 
(the amplitude of the n g = resonance is much greater than that of its "shadow" at n g = 0.5), which suggests that the quantity 
8 A — (1 /2)Eq remains positive (albeit small) at n g = 0. 

The sensitivity of the studied circuit to the charge on the CPT island is illustrated by Fig. 6. Figure 6b shows the time 
dependence of the phase <p of transmitted microwaves when the microwave frequency is tuned to the resonance of the "CPT 
+ resonator" circuit. The noise was measured at full flux frustration (0# = it) when an avoided crossing between the CPT and 
resonator levels was observed, but relatively far from the charge degeneracy point (n g = 0.17). The amplitude of the observed 
telegraph noise corresponds to the charge fluctuations Aq s» 0.03e due to, presumably, coupling of the CPT island to a single 
charge fluctuator in its environment. 



We have performed a detailed analysis of the microwave response of Cooper pair transistors with Ej/Ec ~ 1.5 — 3 coupled to 
a lumped element resonator as a function of the magnetic flux and the gate voltage. Away from the full frustration in flux and 
charge the levels of the Cooper pair transistor are far away from the resonator frequency. In this regime the modulation of the 
resonator frequency induced by CPT can be described as the effective inductance of the CPT that adds to the inductance of the 
resonator. Close to the full frustration the frequency of the resonator approaches and eventually crosses the excitation level of 
the Cooper pair transistor resulting in a complex pattern of the resonator frequency dependence on the flux and gate voltage. In 
all regimes the resonator frequency dependence is very well described by the results of the numerical diagonalization of the full 
Hamiltonian that contains the resonator and CPT. High sensitivity of the resonator frequency close to the level crossing provides 
the tool to measure charge fluctuations in the environment with high accuracy and short time resolution. 
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Figure 5: The color-coded plots of the microwave amplitude 152 il versus the microwave frequency / and the charge n g induced on the CPT 
island by the gate voltage. The phase difference across the CPT was controlled by the magnetic flux <t> in the device loop: (a) $>/$>o = 0.58 
, (b) <1>/<£>o = 0.75, (c) 4>/<&o = 0.90, (d) $>/<t>o = 0.94, (e) "I>/3>o = 1- The solid curves show the numerical simulations with the fitting 
parameters discussed in the text. 
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